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Let —A denote the Dirichlet Laplace operator on a bounded open set in K d . We study 
the sum of the negative eigenvalues of the operator — h 2 A — 1 in the semiclassical limit 
h — > 0+. We give a new proof that yields not only the first term of the asymptotic 
formula but also the second term involving the surface area of the boundary of the set. 
The proof is valid under weak smoothness assumptions on the boundary. 
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1. Introduction and main result 
1.1. Introduction 

Let f2 be a bounded open set in R d , d > 2. We consider the Dirichlet Laplace 
operator —An defined as a self-adjoint operator in L 2 (f2) generated by the form 

(v,-A n v) = / \Vv(x)\ 2 dx 
J ft 

with form domain Hq(Q). Since fl is bounded the embedding of Hq(£1) into L 2 (il) 
is compact and the spectrum of — Aq is discrete. It consists of a series of positive 
eigenvalues < Ai < A2 < . . . accumulating at infinity only. 

In general, the eigenvalues cannot be calculated explicitly and especially 
for large k it is difficult to evaluate them numerically. Therefore it is interesting 
to describe the asymptotic behavior of Afc as k —¥ 00. This is equivalent to the 
asymptotics of the negative eigenvalues of the operator 
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H a = -h 2 A n - 1 

in the semiclassical limit h — > 0+. 

The first general result is due to H. Weyl who studied the counting function 

Nn(h) = U^k<h- 2 } = Tr(Hnt ■ 

In 1912 he showed that the first term of its semiclassical limit is given by the 
phase-space volume [11]: For any open bounded set fi C M. d the limit 

N n (h) = C d \n\h- d + o(h- d ) 

holds as h — > 0+, where 




and uJd denotes the volume of the unit ball in R d . 

H. Weyl conjectured in [12] that this formula can be refined by a second term of 
order h~ d+1 depending on the boundary of fi. This stimulated a detailed analysis of 
the semiclassical limit of partial differential operators. We refer to the books [USE E] 
for general results and an overview over the literature. Eventually, the existence of 
a second term was proved by V. Ivrii by means of a detailed microlocal analysis [5] : 
If the boundary of Vl is smooth and if the measure of all periodic geodesic billiards 
is zero then the limit 

N n (h) = C d \n\ h- d - |0fi| h- d+1 + o(h- d+1 ) (1) 

holds as h — > 0+, where \dCl\ denotes the surface area of the boundary. 
In this article we are interested in the sum of the negative eigenvalues 

Ti-(Hn)- = ]T(/* 2 A fc - 1)_ . 

This quantity describes the energy of non-interacting, fermionic particles trapped 
in n and plays an important role in physical applications. 

The asymptotic relation ([T]) immediately implies a refined formula for the semi- 
classical limit of Tt(Hq)-'- Suppose that the aforementioned geometric conditions 
on f2 are satisfied. Then integrating (JT]) yields 

Tr(Hn)- = L d \Q\ h~ d - ±L d -i\dSl\hr' i+l + o(h- d+1 ) (2) 
as h — > 0+, with 

Ld = I (bl 2 -i)-^= .J,^ 

J Rd d{d + 2) (27r) o! 

In the following we present a direct approach to derive the semiclassical limit of 
Ti(Hq)-. We prove ([2]) without using the result for the counting function. Since we 
do not apply any microlocal methods the proof works under much weaker conditions. 
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1.2. Main Result 

Our main result holds without any global geometric conditions on 51. We only require 
weak smoothness conditions on the boundary - namely that the boundary belongs 
to the class C 1 '™ for some a > 0. That means, we assume that the local charts of 
51 are differentiable and the derivatives are Holder continuous with exponent a. 

Theorem 1.1. Let the boundary of 51 satisfy <9S1 E C 1,a , < a < 1. Then the 
asymptotic limit 



holds as h — > 0+. 

Our work was stimulated by the question whether similar two-term formulae 
hold for non-local, non-smooth operators. This is unknown, since the microlocal 
methods leading to ((T|) are not applicable. Therefore it is necessary to use a direct 
approach. 

Indeed, Theorem 11.11 can be extended to fractional powers of the Dirichlet 
Laplace operator [3] . The strategy of the proof is similar but dealing with non-local 
operators is more difficult and elaborate. In order to give a flavor of our techniques 
we confine ourselves in this article to the local case. 

The question whether the second term of the semiclassical limit of Ti(Hq)- 
exists for Lipschitz domains 51 remains open. 

1.3. Strategy of the proof 

The proof of Theorem 1 1.1 1 is divided into three steps: First, we localize the operator 
Hn into balls, whose size varies depending on the distance to the complement of 51. 
Then we analyze separately the semiclassical limit in the bulk and at the boundary. 

To localize, let d(u) — inf{|a; — u\ : x ^ 51} denote the distance of u G M d to 
the complement of 51. We set 



where < Iq < 1 is a parameter depending only on h. Indeed, we will finally choose 
proportional to h 2 / ( Q + 2 ) . 

In Section [3] we introduce real- valued functions cf> u S C DO (M d ) with support in 
the ball B u — {x £ M. d : \x — u\ < l(u)}. For all u € K d these functions satisfy 




i( u ) = 1 -(i + (d(uf + iiy 1/2 ) 



Halloo < c, 



IIWJL < ci{u) 



(3) 



and for all x E M! 




(4) 
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Here and in the following the letter C denotes various positive constants that might 
depend on J7, but that are independent of u, Iq and h. 



Proposition 1.1. For < Iq < 1 and h > we have 



Tr(H r 



Tv^ u H n ^ u )_l(u)- d du 



< ci^h- d+2 . 



In view of this result, one can analyze the local asymptotics, i. e., the asymptotic 
behavior of Tr((j) u Hn(f) u )- separately on different parts of 51. First, in the bulk, where 
the influence of the boundary is not felt. 

Proposition 1.2. Assume that <\> <E C{f(fl) is supported in a ball of radius I > 
and that 



||V0]|oo < CI- 1 

is satisfied. Then for h > the estimate 

Tr(<f>H n <f>)_-L d [ 4?{x)dxh- d 
Jn 



(5) 



< C 



jd-2 



holds, with a constant depending only on the constant in ([5p. 

Close to the boundary of J7, more precisely, if the support of <f> intersects the 
boundary, a term of order h~ d+1 appears: 

Proposition 1.3. Assume that <\> 6 C§°(M. d ) is supported in a ball of radius I > 
intersecting the boundary of £1 and that inequality {5p is satisfied. 
Then for all < I < 1 and < h < 1 the estimate 



Tr UHc 



L d / (j) 2 (x)dxh d + -L d - 1 I (f) 2 (x)da(x)h 



-d+l 



Oil 



< r(l,h) 



holds. Here da denotes the d— 1-dimensional volume element of dQ and the remain- 
der satisfies 



r(l,h) < C 



yl-2 pa+d-1 j^d+c 



h d-2 h d-l h d 

with a constant depending on f2, Halloo and the constant in {5p. 

Based on these propositions we can complete the proof of the main result. 

Proof of Theorem 11.11 In order to apply Proposition 11.31 to the operators 
(f> u Hn4>u, we need to estimate l(u) uniformly Assume that u € M. d satisfies 
B u n <9f2 ^ 0. Then we have d(u) < l(u), which by definition of l(u) implies 



l(u) < lo/Vz. 



(6) 
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In view of ([3]) we can therefore apply Proposition 11.21 and Proposition 11.31 to all 
functions <f> u , u £ M. d . Combining these results with Proposition ll.il we get 



Tr(# n )_ 



h d 



du 



+ 



l{u) d 4ft d "i J Rd J dn 



4>l(x)da(x) 



du 



l(u) d 



< C [ i^h- d+2 



+ / l{u)- z duh- d+z + / r(l(u),h)l(u)- d du) , 



f/2 



where Ui = {u G Cl : B u ndO, = 0} and U 2 = {u G R d : B u n dQ ^ 0}. Now we 
change the order of integration and by virtue of (jj) we obtain 



Tr (H n )_ - L d \n\ h- d + -L d _x |0fi| 



<C[lo L h- d+2 + l(u)- z duh- d+z + I r(l{u),h)l{u)- d du 



(7) 



It remains to estimate the remainder terms. Note that, by definition of l(u), we 
have 

l( u ) — ^min(d(«),l) and l(u) > 



for all u £M. d . Together with ^ this implies 



l{ U y 2 du < CIq 1 and / l(u) a du < Cl% 



{d(u)<l } 



du < Cl% +1 (8) 



for any a £ R. Inserting these estimates into ([7]) we find that the remainder terms 
are bounded from above by a constant times 

l^h- d+2 + ll a h- d+1 +l% +1 h- d . 

Finally, we choose lo proportional to ft, 2 /( a + 2 ) and conclude that all error terms in 
equal 0(/i- d + 1 + a /( 2 +«)) as h -> 0+. □ 

The remainder of the text is structured as follows. In Section [5] we analyze the 
local asymptotics and outline the proofs of Proposition 11.21 and 11.31 In Section [3l 
we perform the localization and, in particular, prove Proposition ll.il 



2. Local asymptotics 

To prove the propositions we need the following rough estimate, a variant of the 
Berezin-Lieb-Li- Yau inequality (H [5] . 

Lemma 2.1. For any (j> <E C£° (K d ) and h > 

Tr(4>Hn(f>)_ < L d [ cj) 2 {x) dx hr d . 
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Proof. Let us introduce the operator 

H = -h 2 A - 1 , 

denned with form domain ff 1 (R"). The variational principle for sums of eigenvalues 
implies Tr(0_ffo,<^)_ < Tr(0(7?o)_0)_. Using the Fourier-transform one can derive 
an explicit expression for the kernel of (-Hq)_ an d inserting this yields the claim. □ 

2.1. Local asymptotics in the bulk 

First we assume cf> £ Cq°(Q). Then we have Tr (c/>HQ(j>)_ = Tr (cf>Ho(f>)_, since the 
form domains of <^>-ffo,0 and 4>Ho4> coincide. Moreover, by scaling, we can assume 
1 = 1. Thus, to prove Proposition II .21 it suffices to establish the estimate 

Tr(^if ^)_ -L d I cj) 2 {x)dxhr d 

for h > 0. The lower bound follows immediately from Lemma |2. II The upper bound 
can be derived in the same way as in the proof of Lemma 12.31 below. Indeed, by 
choosing the trial density matrix 7 = x(Hq) _x we find 

Tr(0i7 o 0)_ > L d [ <p 2 (x)dx-C d [ {\7<j)) 2 {x)dxh- d+2 
and the claim follows. 

2.2. Straightening the boundary 

Here we transform the operator Hq locally to an operator given on the half-space 
R d ={y£R d :y d > 0}. There we define the operator H + in the same way as i/ji , 
with form domain if^R^). 

Under the conditions of Proposition 11.31 let B denote the open ball of radius 
I > 0, containing the support of <p. Choose xo € B n <9f2 and let v Xo be the normed 
inner normal vector at xq. We choose a Cartesian coordinate system such that 
x = and v Xa = (0, . . . , 0, 1), and we write x = (x', x d ) € x E for x € K d . 

For sufficiently small I > one can introduce new local coordinates near the 
boundary. Let D denote the projection of B on the hyperplane given by x d = 0. 
Since the boundary of £1 is compact and in C ,a , there exists a constant c > 0, such 
that for < I < c we can find a real function / £ C l a given on D, satisfying 

dQ n B = {(x', x d ) : x' £ D, x d = f(x')} n B . 

The choice of coordinates implies /(0) = and V/(0) = 0. Since / £ C 1,Q and the 
boundary of Q is compact we can estimate 

sup |V/(x')| < Cl a , (9) 
x'eD 

with a constant C > depending only on f2, in particular independent of /. 

Now we introduce new local coordinates given by a diffeomorphism ip : Dxl-> 
R d . We set yj — <Pj(x) — Xj for j = 1, . . . , d — 1 and y d — (p d (x) = x d — f(x'). 



< Ch- d + 2 
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Note that the determinant of the Jacobian matrix of tp equals 1 and that the inverse 
of ip is defined on ranip = DxI. There we define </> = <f) ° l P~ 1 and extend it by 
zero to R d , such that $ £ Cq (R d ) and HV^Ik, < Cr 1 holds. 



Lemma 2.2. For < I < c and any h > the estimate 



Tr(0i^ 
holds. Moreover, we have 

and 



Tr(^fT 



(f) 2 (x) dx 



< Cl d+a h- d 



4> 2 (y)dy 



<fi 2 (x) dcr(x) 



4?{y',Q)dy' 



< Cl d ~ 1+2a . 



(10) 



(11) 



(12) 



Proof. The definition of <j> and the fact det Jtp = 1 immediately give (fTTj) . Using 
© we estimate 



4> 2 (y', 0)VT+Wjfdy' < / 4> 2 W, o)dy' + ci d - 1+2a 



4> 2 (x)da(x) = 

on 

from which (fT2")l follows. 

To prove (fTU|) fix v <E -ffp(^) with support in B. For y g rant/3 put = 
u o <^ _1 (y) and extend £ by zero to Kl d . Note that v belongs to Hq^Mtl). 

An explicit calculation shows 

{v, - A K d 5) - («, - A n «) | < C r (5, - A K d 5) . 
Hence, we find 

TitfHnfl- < Tr(0(-(1 - CZ a )/i a A H: , - 1)0)_ . 
Set e = 2Cl a and assume Z to be sufficiently small, so that < e < 1/2 holds. Then 

Tr(0ff o 0)- < Tv(4>(-(l-Cl a )h 2 A K -1)0)- 

< TtC^-^Ar, - 1)£)_ + Tr(0(-(e - CZ Q )/i 2 A R d - e)0)_ 

< Tr(0#+0)_ +eTr(0(-(/i 2 /2)A R d - 1)0)_ . 

By Lemma [Owe have Tr(^(-(/i 2 /2)A R d - 1)0)_ < Cl d h- d and we obtain 

Tr(0ir o 0)_ < Tr(0i7+0)_ +Cl d+a h- d . 

Finally, by interchanging the roles of Hq and H + , we get an analogous upper bound 

□ 



and the proof of Lemma I2T21 is complete 
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2.3. Local asymptotics in half-space 



In view of Lemma 12.21 we can reduce Proposition 11.31 to a statement concerning 
the operator H + , given on the half-space R+. Indeed, to prove Proposition II .31 it 
suffices to establish the following result. 

Lemma 2.3. Assume that <\> G CQ(R d ) is supported in a ball of radius I > and 
that {5p is satisfied. Then for h > the estimate 



Tv^H + ^_-^J^J 2 (x)dx + ^±± I J>\x',Q)dx' 
holds with a constant depending only on the constant in ([5]). 



C Cl d - 2 h- d+2 



Proof. On Ml we can rescale <fi and assume I = 1. In a first step we prove the 
estimate 

Tr^H+cf,) [ <p(x)dx+ [ <t?{x)[ cos(2^ d x d h-')(\e-l)-^- d 

n Jr-^ Jr-^ JR d (znn) 

< Ch- d+2 . (13) 

To derive a lower bound we use the inequality Tr (<fiH + <fi)_ < Tr ((/)(H + )-(j)) 
and diagonalize the operator (_ff + )_, applying the Fourier-transform in the x'- 
coordinates and the sine-transform in the ir^-coordinate. This yields 



Tr(0ff 



+ <P)-<[ ^{x)f 2sin 2 (^0(ie| 2 -l)_ 7 ^ 
Jr± Jr" {2irh) d 



and the lower bound follows from the identity 

2 sin^Xd/i- 1 ) = 1 - cos^x^ 1 ) . (14) 
To prove the upper bound, define the operator 7 = x(-ff + )°_x with kernel 

y) = t^4w x{x) I e l ^*'-y')i h sm{z dXd hr l ) sHUyah-^ x (y) , 

(2nh) d J m<1 

where x denotes the characteristic function of an open ball containing the support 
of (f>. Thus, 7 is a trace-class operator, satisfying < 7 < 1 and by the variational 
principle it follows that 

Tr(0£T + f»_ 
> -Tr(7f>iJ+0) 

,2 \ dC 



= -2 J i (V||Ve l «'-/' l sin(£ d • hr 1 )^^ - \\sm(( d ■ 
>[ (le| 2 -l) _/ ^{x)2s^^ dXd h^)^§- d -Ch-^. 

JR d JR d [ Z7rn ) 

In view of (fT4|) this gives an upper bound and we established (fT3 



\L 2 (Rf) 



(2irh) d 
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We proceed to analyzing the term in (fTB")) which contains the cosine. We substi- 
tute Xd = th and write 

2 (x) / co S (2^ d O(l£| 2 -l)_7^ 

<p 2 (x',th)dx / cos(2^t) (|C| 2 - 1)_ d£,dth- d+1 . (15) 



Note that 



1 Z" 00 f 1 
— / / cos(2£ d i) (|e| 2 - 1)_ de A = jLi-x . (16) 
( 27r ) Jo Jr* 4 

Moreover, in [TJ (9.1.20)] it is shown that 

£ cos(2^) (|C| 2 - 1)_ df = Cjf \oB(2fct)(l - d) (d+1)/2 d& = g ^fgj^ , 

where Jd/2+1 denotes the Bessel function of the first kind. We remark that 
|^d/2+i(2t)| is proportional to t d / 2+1 as t — > 0+ and bounded by a constant times 
t~ 1/2 as £ — y 00 , see [I, (9.1.7) and (9.2.1)]. It follows that 

cos(2&i) (|C| 2 - l)_d£ dt <C t- d / 2 \J d/2+1 (2t)\dt < C. (17) 

Jo 

In view of ([TS]). ([15]) and (JTTJ) we find 

4 * 2(I) L cos(2W "' ) (lf|2 " (W - 5^ *V,<W 

< Ch- d + 2 . 

Inserting this into (fT3f proves Lemma [2~3l □ 
Proposition 11.31 is a consequence of Lemma 12.21 and Lemma 12.31 

3. Localization 

Here we construct the family of localization functions (4> u ) U £R d an d prove Propo- 
sition 11.11 The key idea is to choose the localization depending on the distance to 
the complement of f2, see [4j Theorem 17.1.3] and [TO] . 

Fix a real- valued function <f> g C^(R d ) with support in {|x| < 1} and ||0||2 = 1. 
For u,x £ M. d let J(x,u) be the Jacobian of the map u n- (x — u)/l(u). We define 



4>u{x) =4>[j^j ^/J(x7u)l(u) d/2 , 



such that <p u is supported in {a; : \x — u\ < l(u)}. According to |10j . the functions 
4> u satisfy © and Q for all u G M. d . 

To prove the upper bound in Proposition II .![ put 

7 = / 4> u (4> u Ha(j) u f_ 4> u l{u)~ d du . 
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Obviously, 7 > holds and in view of @ also 7 < 1. The range of 7 belongs to 
Hq(CI) and by the variational principle it follows that 

-Tr(H n )_ < TrjHn = - f Tr (<t> u H n cj> u ) _ l(u)~ d du . 

JR d 

To prove the lower bound we make use of the IMS-formula 

\ (/, 2 (-A)/) + \ (/, -A0 2 /) = (/, 0(-A)0/) - (/, /(V^) 2 ) , 
valid for <p £ C£°(R d ) and / e Hq(Q). Combining this identity with g} yields 

(f,-Af)= [ ((f, ( j )u (~A)4> u f)-(f,^cp u ) 2 f))l{uy d du. (18) 

JR d 

Using ((3]) and (|4]) one can show [10] 

f (V(pu) 2 (x)l(u)- d du < C I 4>l(x) l{u)~ d - 2 du . 

JR d JR d 

We insert this into (jl8]l and deduce 

Tr(H Q )_ < [ r TY{4> u {-h 2 A-l-Ch 2 l{u)- 2 )4>u)_l{u)- d du, 
Jo,* 

where fl* — {u g R d : supp</> u n f2 ^ 0}. To estimate the localization error we use 
Lemma l2.ll For any Hgl, let p u be another parameter < p M < 1/2 and estimate 

Tr (cj> u (-h 2 A - 1 - Ch 2 l{u)- 2 ) cp u )^ 

< Tr (M-h 2 A - 1)0„)_ + CTr (- Att ^ 2 A - p u - /^(u)" 2 ) 4> u ) _ 

< Tr {<f> u Hn<l>u)- + C l{u) d ( fiu h 2 )- d l 2 (p u + h 2 l(u)- 2 ) 1+d/2 . 
With p u proportional to h 2 l(u)~ 2 we find 

Tt(H Q )_ < [ r Tr{<}> u Hn(t>u)^l{u)- d du + Ch- d+2 [ l(u)- 2 du. 

In view of the last integral is bounded by a constant times Iq 1 and the proof of 
Proposition |TTT| is complete. 
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